OCEAN SURFACE CURRENT MAPPING WITH BISTATIC HF RADAR 



Reference to Related Applications 
5 This application claims priority to U.S. Provisional Applications Serial Nos. 

60/308,235 filed July 26, 2001, 60/315,567 filed August 28, 2001, and 60/327,217 filed 
October 4, 2001, all of which are incorporated herein by reference. 

Field of the Invention 

10 The present invention relates generally to radar systems and more particularly to 

a radar systems and method for mapping surface currents in large bodies of water. 

Back ground of the Invention 

Low-frequency backscatter radar systems, operating in the MF, HF, VHF, and low 
15 UHF bands are widely used for mapping surface water velocity such as currents on the 
oceans and flow along rivers. 

The most commonly used are backscatter radar systems that use Doppler 
principles to measure current velocity. Two effects contribute to the backscatter Doppler 
shift: (i) the motions of waves that Bragg scatter, i.e., whose lengths are exactly half the 
20 radar wavelength; and (ii) the component of the underlying current that transports the 
surface waves along the bearing of the radar. For a backscatter radar, the wave Doppler 
occurs at a single frequency and is precisely known, and can be readily subtracted from 
a measured Doppler shift to get the desired radial current speed. Thus, the waves serve 
as tracers of the underlying currents. That is, information about these waves is not being 
25 sought; it is only necessary that they be present to serve as radar targets. 
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Range or distance to the scattering cell is obtained from the time delay between 
transmit and received radar signals as described, for example, in commonly assigned U.S. 
Patent No. 5,361,072, which is incorporated herein by reference. Following range 
processing, the echo time series for each range cell is Fourier transformed to obtain 
5 Doppler spectra and cross spectra among several receive antennas or elements. The Bragg 
peaks used for current extraction are identified and isolated. At each Doppler spectral bin 
(corresponding to a definable current radial velocity), the angular bearing of the echo is 
obtained from the directional properties of the receive antennas using a bearing 
deterrnining algorithm One suitable bearing determining algorithm, for example, is the 
10 Multiple Signal Classification (MUSIC) direction-finding algorithm as described in 
commonly assigned U.S. Patent No. 5,990,834, whichis incorporated herein by reference. 
Thus, the backscatter radar makes its measurements in a polar coordinate system in which 
radial current velocity is mapped as a function of distance and bearing angle. 

A single backscatter radar system can measure only the radial component of a 
15 two-dimensionalhorizontal current velocity vector. Therefore, normally two backscatter 
radar systems are used in pairs, spaced tens of kilometers apart along the coast and 
operating independently. Based on the known geometry and location of a mutually 
observed scattering cell, two resulting radial velocity components are combined to 
produce a total velocity vector map across the overlapping coverage zone. Thus, one 
shortcoming of conventional systems for mapping surface water velocities is the need for 
multiple, costly backscatter radar systems. 

Another shortcoming with conventional systems is that backscatter radar systems 
cannot resolve a total vector on or near a line joining the two backscatter radar systems, 
because both are measuring the same component of velocity. This is particularly a 
25 problem for regions falling along the coast or across rivers and mouths of bays, which are 
often of great interest. 

An alternative approach for mapping surface water velocities uses a bistatic radar 
system in which one or more receivers operate with two or more transmitters spaced apart 
therefrom Resultant echo sets based on the signals simultaneously radiated from each 
30 isolated transmitter-receiver pair, can then be used to construct non-parallel components 
of velocity at the scattering cell. A major drawback to this approach is the expense of 



20 
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synchronizing a geographically separated transmitter and receiver. Previous bistatic 
systems have employed highly stable Cesium or Rubidium timing standards, or 
oven-controlled phase-lock loops to maintain coherent signal bases at the separated 
points. In addition, to being expensive such methods are ill suited to harsh environments 
5 to which ocean sensing radar systems are routinely exposed. 

A more fundamental problem with the use of bistatic radar systems is that in 
contrast to backscatter radar systems, contours of constant time delay for bistatic radar 
systems define ellipses confocal about the transmitter and receiver pair. As a result, the 
Doppler shift due to the wave motion is no longer a constant as it is with backscatter radar 
10 systems, but varies with position around the cell or contours of constant time delay. 

Accordingly, there is a need for a radar system and method of operating the same 
that improves coverage area and accuracy in regions where geometry of the system would 
normally preclude monitoring with conventional backscatter radar systems. It is also 
desirable that the radar system and method be robust and inexpensive. It is also desirable 
1 5 that the radar system and method enable existing ocean sensing radar systems to be easily 
and inexpensively upgraded. 

The radar system and method of the present invention provides these and other 
advantages over the prior art. 

20 .Summary of the invention 

The present invention provides a bistatic radar system and method for mapping 
surface currents in bodies of water that offers both expanded coverage area and stable 
estimates of current velocity in zones or regions where conventional backscatter radar 
systems have not been effective. 
25 It is an object of the present invention to provide oceanic current information using 

a bistatic radar system that improves coverage area and accuracy in regions where 
instabilities normally preclude monitoring with conventional backscatter radar systems. 

It is a further object of the present invention to reduce complexity, cost, and use 
of radio spectral resources by synchronizing the modulating signals of transmitters and 
30 receivers of a bistatic radar system with Global Positioning System (GPS) timing signals. 

It is a yet further object of the present invention to provide a method and computer 
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prograi-n product for transforming or converting time-delay or range, bearing, and velocity 
relationships used in conventional backscatter radar systems into the required current 
velocity maps used with elliptic/hyperbolic geometries of bistatic radar systems. 

According to one aspect of the present invention, a bistatic radar system is 
5 provided having a number of transmitters and receivers for transmitting and receiving 
radar signals. Preferably, at least one transmitter is positioned in a location separate from 
at least one receiver, and each of the transmitters and receivers include a local oscillator 
locked to Global Positioning System (GPS) signals to provide the necessary coherency 
between the transmitters and receivers. More preferably, the transmitters and receivers 
10 are configured to detect and measure oceanic conditions, and the bistatic radar system 
further includes signal processing means adapted to derive information on the oceanic 
conditions. Oceanic conditions detected and measured by the number of transmitters and 
receivers can include, for example, surface current velocity vectors. 

Generally, the signal processing means is adapted to determine a current velocity 
15 within a scattering cell using a Doppler shift (f D ) of the radar echoes, and the following 
equation: 

Jo ±>fe cos f 

where X is the wavelength of the radar signals; g is the acceleration of gravity 
(9.806 m/s 2 ); 0 is a bistatic angle between lines connecting the transmitter and the 

20 scattering cell and the receiver and the scattering cell; and V h is the velocity component 
measurable from the bistatic radar system The current velocity, V* is determined along 
a hyperbola perpendicular to an ellipse passing through the scattering cell and confocal 
about the transmitter and receiver, the ellipse having a constant time delay (D) equal to 
the measured radar echo time delay. In one embodiment, the Doppler shift (f D ) is 

25 measured directly using the bistatic radar system In another embodiment, the signal 
processing means is adapted to calculate a Doppler shift (f D ) using the following equation: 



V 7Li X 



where X is the wavelength of the radar signals; g is the acceleration of gravity 
A-70836/AJTAVEN December 19 , 2001 (1035432) 



(9.806 m/s 2 ); and V r is a pseudo-radial current velocity derived in a pseudo backscatter 
stage using a computer program developed for a backscatter radar system and oceanic 
conditions detected and measured by the plurality of transmitters and receivers. 

Preferably, the bistatic radar system is adapted to provide surface current velocity 
5 vectors that are independent of motions of waves having velocities over a Doppler 
spectral region substantially the same as the surface current velocity vectors. More 
preferably, the bistatic radar system is adapted to provide total current vectors in regions 
along lines joining multiple receivers. 

According to another aspect of the present invention, a method is provided for 
10 mapping surface current vectors using a radar system including a number of transmitters 
and receivers, including a bistatic radar subsystem having at least one transmitter 
positioned in a location separated from at least one receiver. Generally, the method 
involves steps of: (i) scattering radar signals from the transmitter off ocean waves in a 
scattering cell to produce echoes in the receiver; (ii) deteraiining a bearing angle ((p) to 
15 the scattering cell using using a bearing determining algorithm; (iii) sampling versus time 
after transmission to measure radar echo time delay from the transmitter to the receiver; 
(iv) deterrnining the location of the scattering cell; and (v) detenriining a current velocity 
at the scattering cell. 

Generally, the step of determining a current velocity at the scattering cell is based 
20 on the bistatic angle between lines connecting the transmitter and the scattering cell and 
the receiver, the Doppler shift (f D ) of the radar echoes and the position of the scattering 
cell. As provided above, the Doppler shift can either be measured directly from the 
Doppler shift of a bistatic transmitter - receiver pair, or calculated from a hypothetical 
radial backscatter velocity component measured relative to one of the receivers. 
25 The procedure for detenriining a position of the scattering cell includes the steps 

of: (i) determining a major axis (A) of an ellipse passing through the scattering cell and 
confocal about the transmitter and receiver, the ellipse having a constant time delay (D) 
equal to the measured radar echo time delay; (ii) determining a minor axis (B) of the 
ellipse; (iii) deterrnining sine and cosine of an angle (Y) from the scattering cell to an 
30 origin of a local coordinate system; and (iv) detennining the position of the scattering cell 
within the local coordinate system from the sine and cosine of the angle and the major and 
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minor axis of the ellipse. 

In yet another aspect, the present invention is directed to a computer program 
product for use in conjunction with a computer system to accomplish steps of the above 
method. The computer program product includes a computer readable storage medium 

5 and a computer program mechanism embedded therein. The computer program 
mechanism, includes a program module that directs the computer system, to function in 
a specified manner, to map surface current vectors using a radar system including a 
number of transmitters and receivers, including a bistatic radar subsystem having at least 
one transmitter positioned in a location separated from at least one receiver. Generally, 

10 the computer program includes: (i) a radar control subroutine or program module; (ii) a 
bearing angle (cp) detenriining subroutine or program module; (iii) an echo time delay 
subroutine or program module; (iv) a position deter m ining subroutine or program module; 
and (v) a current velocity subroutine or program module. 

Advantages of the computer program and method of the present invention include 

15 any one or all of the following: 

(i) reduced complexity, cost, and use of radio spectral resources through 
synchronization of transmitters and receivers using GPS timing signals; 

(ii) ability to quickly and relatively inexpensively improve the coverage area 
and accuracy of conventional current-mapping backscatter radar systems 

20 by deploying additional transmitters and/or receivers; 

(iii) compact transmit subsystem designs, requiring no computer nor air- 
conditioning resources, and suitable for remote deployment where solar 
energy is a preferred source; 

(iv) ability to radiate from vertical whip transmit antennas that can be mounted 
25 on buoys, offshore structures, or building rooftops without incurring the 

resulting complications of antenna pattern distortions that can affect the 
ability to measure bearing angles accurately; and 

(v) elimination of antenna pattern distortions caused by uncontrolled rotation 
of an antenna mounted on a buoy or offshore structure due to wave 

30 motion. 
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Brief Description of the Drawings 

The invention will be more clearly understood from the following description 
when read in conjunction with the accompanying drawings of which: 

FIG. 1 is a block diagram showing the hardware components and signal processing 
5 stages of current information from bistatic radar systems, according to an embodiment of 
the present invention; 

PIG. 2 is a plan view showing the location of transmitter and receivers of a bistatic 
radar system according to an embodiment of the present invention in and around a body 
of water; 

10 FIG. 3 is a graph illustrating relationships between the measurable current vectors 

and the geometry of both backscatter and bistatic radar systems according to an 

embodiment of the present invention; 

FIG. 4 is a flowchart showing a method of operating a bistatic radar to map surface 

current vectors according to an embodiment of the present invention; 
15 FIG. 5 is a flowchart showing a method for determining the position of the 

scattering cell according to an embodiment of the present invention; and 

FIG. 6 is a block diagram of an embodiment of a computer program showing 

modules for mapping surface current vectors according to an embodiment of the present 

invention. 

20 

Detailed Description 

The present invention is bistatic radar system and a method for operating the same 
to map surface currents in bodies of water, such as oceans, seas, bays, harbors, rivers and 
lakes. 

25 A bistatic radar system according to an embodiment of the present invention will 

now be described with reference to FIG. 1. FIG. 1 is an exemplary block diagram of a 
bistatic radar system 100 including a number of transmitters and receivers for transmitting 
and receiving radar signals. At least one transmitter is positioned in a location separated 
from at least one receiver, and the necessary coherency between the transmitters and 

30 receivers provided using a Global Positioning System (GPS) signal according to an 
embodiment of the present invention. For purposes of clarity, many details of radar 
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systems that are widely known and are not relevant to the present invention have been 
omitted. 

FIG. 1 shows the hardware components and signal processing stages of the bistatic 
radar system 100. It is noted that blocks shown as sharp-cornered rectangles represent 

5 hardware components. Rounded-corner blocks represent the key signal processing steps 
performed on the received radar data, and can be executed by a software program running 
in real time on a digital computer. 

Referring to FIG. 1, the bistatic radar system 100, generally includes a transmitter 
102 that cooperates with a receive- section 104 having several individual receivers 106 to 

10 enable investigation of oceanic or water surface conditions. Receive- antennas 108 of the 
receivers 106 receive echoes of radar signals transmitted from a transmit- antenna 1 10 of 
the transmitter 102 and scattered from waves in a scattering point or cell (not shown in 
this figure). A signal processor 1 12 processes information derived from the radar echoes 
to deterrnine and display information on oceanic conditions. Such information includes 

15 surface current velocity vectors. To operate in the bistatic mode the transmit-antenna 1 10 
is geographically separated from at least one receive-antennas 108, as shown in FIG. 2. 
That is, the transmit-antenna 110 and at least one receive-antenna 108 are located so that 
an angle between the line joining the transmit-antenna and a scattering cell 114, and the 
line joining the receive-antenna and the scattering cell is non-zero. Preferably, to 

20 maximize the area covered by the bistatic radar system 100, all of the receivers 106 or 
receive-antennas 108 are separated from the transmit-antenna 110 and from one another. 
No physical connections or wires link or provide a communication path between the 
transmitter 102 and receivers 106. Distances between the two depend on the frequency 
and mapping coverage needed from the bistatic radar system 100 and can range from 

25 about 2 to about 200 kilometers (km). In a backscatter radar system 1 1 6, for contrast and 
illustration, the transmitter and transmit-antenna are co-located at the same position as the 
receiver. 

It will be appreciated that, the electronic components of the transmitter 102 and 
the receivers 106 need not be near their respective antennas 108, 1 10. Nor do receivers 
30 106 need to be near to the signal processor 1 12. For example, it is often advantageous to 
place the signal processor 112, from which a user or operator accesses information on 
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ocean conditions, in a shelter some distance from the receivers 106. The signal processor 
1 12 can be linked to the receivers 106 via a land-line or a wireless communication path. 
Similarly, the electronic components of the transmitter 102 and the receivers 106 can be 
in an enclosure or building to which the antennasl08, 1 10 are affixed, and coupled to the 
5 antennas via a short cable or wire. 

Advantageously, the transmitter 102 and the transmit-antenna 110 are a self- 
contained unit powered, for example, using photovoltaic cells, a fuel cell, or battery. In 
one preferred embodiment, the transmitter 102 is placed on a water-borne platform such 
a boat, barge, platform or buoy, and surrounded partly or entirely by a ring of receive- 
10 antennas 108. For example, the transmitter 102 can be on a buoy in a harbor or bay 1 1 8, 
and the receive- antennas 108 can be on other buoys or shore installations such as 
buildings on land 120 around the bay, as shown in FIG. 2. 

Generally, the transmitter 102 includes a digital synthesizer 122 having a local 
oscillator (not shown) for generating a carrier frequency or wave and a signal for 
15 modulating the carrier wave to produce the radar signal transmitted by the transmitter. 
For a bistatic configuration the radar signal can take the format of, for example, a 
continuous 20 MHz carrier wave modulated as described in U.S. Patent 5,361,072. The 
modulated radar signal is then passed from the digital synthesizer 122 through a transmit- 
amplifier or power-amplifier 124 and radiated from the transmit-antenna 110. The 
20 transmit-antenna 110 is usually configured to provide a broad beam, usually between 
120°and 360° in bearing angle, illuminating a large area of ocean or water concerning 
which information on oceanic conditions is needed. 

The receivers 106 include a receiver module 126 for receiving and decoding radar 
echoes, and several signal processing stages or signal processors for extracting 
25 information the received echoes. The receiver module 126 mixes the received echoes 
with a signal generated by a local oscillator (not shown) to decode the echoes. Generally, 
the signal processors have both hardware and software elements. Preferably, each 
receiver module 126 includes an analog-to-digital (A/D) convenor (not shown) to provide 
an output that is a digital time series that can be digitally processed by the signal 
30 processors. A digital signal processor 128 determines the position of a time-delay cell and 
outputs a time series of radar echo and noise for each range cell. A Doppler signal 
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processor 130, hi a process known as coherent integration or processing, orders potential 
target echoes into a two-dimensional space of range and Doppler frequency. A final 
processing stage, a bearing angle processor 132, identifies signals of interest and separates 
them from background noise and clutter before deteraiining an echo bearing angle. Each 
of these three signal processors 128, 130, 132, are described in greater detail in U.S. 
Patents 5,361,072 and 5,990,834, which are incorporated herein by reference. 

Because precise timing is required between the radiated signal and the received 
signals to allow coherent radar processing, the transmitter 102 and each receiver 106 
further include a GPS satellite receiver 134 to synchronize the modulating signal in the 
transmitter with the locally generated signal in the receiver module 126. These GPS 
satellite receivers 134 are small devices that come equipped with their own satellite 
antenna 136 and are commercially available from several vendors, including, for example, 
the Thuderbolt GPS Disciplined Clock sold by Trimble Navigation Limited of Sunnyvale, 
California. One method of utilizing a GPS timing signal for synchronizing the modulating 
signals of the transmitter 102 and receivers 106 is described, for example, in commonly 
assigned, co-pending U.S. Provisional Applications Serial No. 60/315,567 filed August 
28, 2001, and incorporated herein by reference. Local phase locked oscillators (PLO) in 
the transmitter 102 and receivers 106 lock to a GPS signal from the GPS satellite 
receivers 134 to provide a common reference. The PLO serves as an electronic flywheel 
that over a period of time builds up to an accuracy of 10 12 to 10 13 parts per million. 

The signal processor 1 12 includes a general purpose digital computer, such as a 
personal computer or PC, programmed according to the present invention to process 
information derived from the radar echoes to determine and display information on 
oceanic conditions. Generally, the signal processor 112 encompasses several elements 
or stages in which information is processed, including a bearing angle determination stage 
138, a bistatic scattering point determination stage 140, a bistatic current velocity stage 
142, and a total current velocity vector stage 144. Eachof these stages and the operations 
performed therein will now be described with reference to FIGs. 1 and 3. 

In the first or bearing angle determination stage 138 the bearing or bearings, range 
and Doppler for each scattering cell 1 14 or cells are determined from digital target echo 
signals received from each receiver 106. In a preferred embodiment, this determination 
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made using a suitable bearing determiniiig algorithm, such as the Multiple Signal 
Classification (MUSIC) direction-finding algorithm described in U.S. Patent No. 
5,990,834. This stage, the bearing angle determination stage 138, can also be called a. 
pseudo-backscatter stage, since, as in a backscatter radar system in which a transmitter 
5 and receiver co-located, the output from this stage 138 can be represented as a polar 
mapping of radial surface current velocity versus range and bearing. Each Doppler shift 
(or spectral bin) is related to the velocity component of the scattering cell 114 moving 
toward or away from the backscatter radar system This by definition is the radial velocity 

m. 

10 In one embodiment, the output of thepseudo backscatter stage 138, is transformed 

into a bistatic total current velocity vector map corresponding to the radial surface current 
velocity map of a backscatter radar system The mathematics and physics of this 
transformation, from a radial surface current velocity map to a total current velocity vector 
map, will now be described with reference to FIG. 3. Referring to FIG. 3, the transmitter 
15 102 of a bistatic radar system 100 according to the present invention is located at a circle 
146 at a position x,y= 14.1, 14.1 of a coordinate system The location of the receiver 106 
is shown by circle 148 at the origin, x, y = 0, 0, of the coordinate system For purposes 
of illustration and comparison this is also the location of a monostatic or backscatter 
transmitter that forms a backscatter radar system when operating with this receiver 106. 
20 For the backscatter radar system, contours of constant time delay between transmission 
and reception of the echo signal to targets form circles 150, 152, concentric to the location 
of the backscatter radar system, shown here as circle 148. In contrast, for the bistatic 
radar system 100, contours of constant time delay define a family of ellipses 154, 156, 
confocal about the transmitter 102 at circle 146 and receiver 106 at circle 148. 
25 Consequently, the interfocal distance of the family of ellipses, denoted here as F, the 
distance between transmitter 102 and receiver 106, and in FIG. 3 is equal to twenty units. 

There are three output observables for backscatter radar systems, whose analogue 
bistatic counterparts are sought to define the desired bistatic total current velocity vector 
map: (i) range to the scattering cell 114; (ii) bearing angle, <p, to the scattering cell; and 
30 (iii) the radial velocity, V r of the water (the current) at the scattering cell. The 
corresponding outputs or counterparts for the bistatic radar systems include: (i) position 
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of the scattering ceU 1 14 (two coordinates are required to define horizontal position); and 
(ii) the current velocity component available from the bistatic radar at the scattering cell. 

A local coordinate system (not shown) is employed for the bistatic configuration, 
whose origin is the midpoint of the line along F between the receiver 106 and transmitter 
102. The x' axis of this local coordinate system is directed toward the transmitter 102 at 
circle 146, and the y" axis is counterclockwise perpendicular to it. 

The position of the scattering cell 1 14 at R, x\ y\ in the local coordinate system 
is determined in the bistatic scattering point determination stage 140 using the following 
mathematical steps. 

The measured radar echo time delay, D, of an ellipse passing through the 
scattering cell 1 14 is related to the distance from receiver 106 to scattering cell, P, the 
distance from scattering cell to the transmitter 102, Q, and the distance, F, separating the 
transmitter and receiver 106 by the following equation: 

D=P + Q-F (Eq. 1) 

The major axis, A, of the ellipse passing through the scattering cell 1 14 with the 
given time delay D is given by the following equation: 

A=D + F CBq.2) 
The minor axis, B, of the ellipse passing through the scattering cell 1 14 with the 
given time delay D, lying along the y" axis is given by the following equation: 

B =^D 2 +2DF (Eq- 3) 

Recalling that the bearing angle, (p, to the scattering cell 114 was measured and 
is known from the bearing angle determination stage 138 using the direction-finding 
algorithm, the following equation provides a denominator, Den, needed for subsequent 
equations: 

Den =(B sec <p) 2 +(F tan cpf (Eq. 4) 

The cosine of an angle from the origin of the local coordinate system, x' = 0, y' = 

0, to the scattering cell 114 at R, is given the following equation: 

— F{F + D) tan* q ) + B 1 sec<? (Eq. 5) 

cosy/ = 

Den 

The sine of an angle from the origin of the local coordinate system, x" = 0, y' = 0, 
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to the scattering cell 1 14 at R, is given the following equation: 



sin y/ = 



(F +(F +D)stc<p)B tan<p (Eq 6) 



Desired position of scattering cell 1 14 in the local coordinate system are given by 
the following equations: 

x'= — cosyr, )>' = — siny/; (Eq- 7) 

2 2 

Current velocity (VJ is along a hyperbola perpendicular to an ellipse passing 
through the scattering cell is determined in the bistatic current velocity stage 142 using 
the following mathematical steps. 

Abackscatter Doppler shift/ D of the echo signal seenby a receiver-transmitter pair 
at circle 148 in FIG. 3 is related to a wavelength ofthe radar and the radial component of 
the current velocity. It is noted that the Doppler shift,/ D , was previously measured in the 
Doppler signal processor 1 30 as described above. The Doppler shift/ D of the echo signal 
from the backscatter configuration is expressed by the following equation: 



fo=± 



_ +2 ^ (Eq-8) 

7iX A 



15 where A. is the radar wavelength, which is known; g is the acceleration of gravity (9.806 
m/s 2 ) and V r , is the radial component ofthe current velocity, directed toward or away from 
the radar. The ± sign denotes the two spectral Bragg peaks seen by HF radars, due to 
waves half the radar wavelength moving toward (+) and away from(-) the radar. Both are 
used in extracting current information, providing redundant estimates, which improves 

20 the accuracy. 

A bistatic Doppler shift f D ofthe echo signal seenby a receiver 106 at circle 148 
and transmitter 102 at circle 146 in a bistatic configuration is expressed by the following 
equation: 

/d=± JT^ +2 ^cos! (Eq.9) 
Jd \nX 2 A 2 

25 where X is the wavelength of the radar signals; g is the acceleration of gravity (9.806 
m/s 2 ); # isabistatic angle between lines connecting the transmitter 102 and the scattering 

A-70836/AJTAVEN December 19, 2001 ( 1035432) 



- 14- 

cell 114, and the receiver 106 and the scattering cell; and V h , is the current velocity 
measurable by the bistatic radar system 100. 

As with the backscatter Doppler shift, Eq. 8, the first term in Eq. 9 represents the 
Doppler shift due to the wave motion. That is, the scatter from Bragg waves advancing 
5 toward or moving away from a baseline separating the transmitter 102 and receiver 106, 
whichhave increasing or decreasing constant time delays, D. This term also suggests the 
complication of bistatic configurations. The Doppler shift due to the wave motion is no 
longer a constant as for backscatter configurations, but varies with position around the 
range cell or ellipse of constant time delay. Thus, it is necessary to separate the Doppler 
10 shift due to the wave motion from the Doppler shift due to current. It is the difficulty of 
this separation that has rendered use of bistatic radar systems to monitor oceanic 
conditions non-obvious in the prior art. The velocity component measurable from the 
bistatic system is called V h in contrast with the backscatter geometry, because it lies along 
a hyperbola that is perpendicular to the ellipse at the scattering point. This vector, V h , is 
15 shownatpointRinFIG.3. As a check of the validity of Eq. 9, note that as the transmitter 
102 and receiver 106 move together so that the configuration becomes that of a 
monostatic or backscatter radar, that is as the bistatic angle goes to zero , Eq. 9 reduces to 
Eq. 8. 

Current velocity (V0 along the hyperbola will now be derived from the measured 
20 Doppler shift using the following equation: 



fp —A/4. 



COSi 



icosf 



(Eq. 10) 



where X is the wavelength of the radar signals; g is the acceleration of gravity (9.806 
m/s 2 ); T? is a bistatic angle; and V h , is the current velocity along the hyperbola passing 
through the scattering cell 1 14. 

Next, the direction of V h with respect to the local x 1 , y 1 coordinate system is 
determined as follows. Recalling that the bistatic Doppler shift was previously measured 
in the Doppler signal processor 130, the steps below show how to determine the required 
bistatic angle. 

The following equations define vectors P, Q and their unit vectors, p' , q' that point 
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away from the receiver and transmitter, respectively toward the scattering point: 

P = {.V + F 12, ?}, p'= j£j Q = {x'-F 12,?}, <f = j|j (Eq. 11) 

Normal and unit normal vectors outward from the ellipse at the scattering point are 
defined by the following equations: 

5 



cos<^ysin<pj; ft=^ (Eq. 12) 



10 



From these sets of unit normal vectors, the cos of (tf /2), which is needed to 
determine the velocity component in Eq. 9, is found using the following equation: 

cos(f) = p'-n = *.n (Eq-13) 

Substituting the cos of ( V IT), into Eq. 9 provides the magnitude of bistatic vector 
V h , and its components within the local coordinate system along x\ y' -- defined as u', v*~ 
are now given as: 

(u',v')=v h -n (Eq-14) 

These are the final current velocity components measurable by the bistatic radar 
system 100, at the point x', y' in the local coordinate system defined by the backscatter or 
transmitter-receiver geometry. Based on the many Doppler spectral bins, a map of 
hyperbolic velocity components, V h , can now be constructed from the measured sea echo 
data. 

It will be appreciated that these velocity components u\ v' at points x\ y' defining 
) this transmitter-receiver oriented map can be transformed and rotated to any other 
coordinate system that is convenient, by well known mathematical techniques. For 
example, if the coastline is oriented along the heavy dashed line 158 of FIG. 3, then the 
coastline-oriented system shown might be more appropriate. Or in other cases, an 
earth-based latitude-longitude (East-North) system is selected. 
5 The hyperbolic velocity, V h , based current mapping derived from the bistatic radar 

system 100 as outlined above is not a complete representation of horizontal current flow. 
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Neither is the radial current map derivable from a backscatter radar system having a single 
backscatter radar. Both systems measure only one component of the two-dimensional 
current vector at each point on the map. Therefore a total current velocity vector map is 
determined in a final stage, the total current velocity vector stage 144. 
5 Prior art backscatter radar systems required at least two complete backscatter 

radars viewing the same point on the surface from two different directions to produce a 
radial current vector map. Moreover, as noted above, this precluded determining current 
vectors at and near the "baseline" joining the two backscatter radar systems. For example, 
if a second backscatter radar were located along the heavy dashed line 158 representing 
10 the coast in FIG. 3, say at position (20,0), the backscatter radar system could not produce 
current vector maps covering the important area near the coast where the angle between 
the two vectors approaches zero. 

In contrast, the bistatic radar system 100 of the present invention having a 
receiver-transmitter pair at circle 148 augmented by a transmitter 102 at circle 146, can 
1 5 produce an accurate total current vector estimates over the full area covered by the bistatic 
radar system, including the area along the coast. Circles define the constant time delay 
contours (or range cells) at which radial maps from the receiver-transmitter pair at circle 
148 operating in a backscatter form produces radial vectors, V r . Operating 
simultaneously, measurements from the bistatic configuration of transmitter 102 at the 
20 circle 1 46 and a receiver 106 at the circle 148 produce hyperbolic vectors, V h , at the same 
coastal points. As shown, these vectors are no longer parallel. In fact, they lie at a 
considerable angle to each other. Thus, a bistatic radar system 100 according to the 
present invention, having two transmitters 102, a single receiver 106, and a signal 
processor 112, is capable of providing a total current velocity vector map over the full 
25 area covered by the system, including the area along the coast. 

A method or process for mapping surface current vectors according to an 
embodiment of the present invention will no w be described with reference to FIG. 4. FIG. 
4 is a flowchart showing steps of a method for operating a bistatic radar system according 
to an embodiment of the present invention. In the method, radar signals are radiated from 
30 the transmitter and are reflected off a scattering cell to produce echoes in the receiver (step 
160). A bearing angle (<p) to the scattering cell is determined using a bearing determining 
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algorithm, such as the MUSIC direction-finding algorithm or any other suitable bearing 
angle determining algorithm (step 162). The radar echoes are sampled versus time after 
transmission to measure radar echo time delay from the transmitter to the receiver (step 
164). The position of the scattering cell is determined (step 166), and the velocity of the 

5 surface current at the scattering cell determined or calculated using Eq. 10 given above 
(step 168). Generally, the step of determining current velocity at the scattering cell, step 
168, involves determining a Doppler shift (f D ) for the radar echoes. In one embodiment, 
this Doppler shift is measured directly using the separated transmitter and receiver of the 
bistatic radar system Alternatively, the Doppler shift may also be calculated from the 

10 bearing angle (<p) and a radial current velocity relative to the receiver, V r , using Eq. 8 
given above. 

A method or process for deteraiining the position of the scattering cell according 
to an embodiment of the present invention will now be described with reference to FIG. 
5. Referring to FIG. 5, the method can be accomplished by: (i) deteraiining a major axis 
1 5 (A) of an ellipse passing through the scattering cell and confocal about the transmitter and 
receiver, the ellipse having a constant time delay (D) equal to the measured radar echo 
time delay (step 170); (ii) detennining a minor axis (B) of the ellipse (step 172); (in) 
deterrnining sine and cosine of an angle (Y) from the scattering cell to an origin of a local 
coordinate system (step 174); and (iv) determining the position of the scattering cell 
20 within the local coordinate system from the sine and cosine of the angle and the major and 
minor axises of the ellipse (step 176). Generally, the step of determining the major axis, 
step 170, is accomplished using Eq. 2 given above. The step of determining the minor 
axis, step 172, is accomplished using Eq. 3, and the step of determining sine and cosine 
of Y, step 174, is accomplished using Eqs. 4 and 5. Finally, the step of deterrnining the 
25 position of the scattering cell of the within the local coordinate system, step 176, is 
accomplished using Eqs. 7 above. 

A computer program for performing each step of the above processes will now be 
described with reference to FIG. 6. FIG. 6 is an illustrative block diagram of a computer 
program 178 according to an embodiment of the present invention. The computer 
30 program 178 generally includes: (i) a radar control subroutine or program module 180; 
(ii) a bearing angle ((p) determining subroutine or program module 182; (iii) an echo time 
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delay subroutine or program module 184; (iv) a position detennining subroutine or 
program module 186; and (v) a current velocity subroutine or program module 188. The 
radar scattering program module 1 80 includes program code adapted to operate or control 
the bistatic radar system 100 to scatter radar signals from the transmitter 102 off the 
5 scattering cell 114 to produce echoes in a receiver 106. The bearing angle ((p) 
determining module 1 82 includes program code adapted to determine a bearing angle ((p) 
to the scattering cell 1 14 using a bearing detennining algorithm The echo time delay 
program module 184 includes program code adapted to measure a time delay of the 
echoes. The position detennining module 186 includes program code adapted to 
10 determine the position of the scattering cell 1 14 within the local coordinate system The 
current velocity program module 188 includes program code adapted to detennining the 
cunent velocity at the scattering cell 1 14. 

In one embodiment, the cunent velocity program module 188 includes program 
code using a Doppler shift measured directly using a geographically separated transmitter 
15 1 02 and receiver 1 06 of the bistatic radar system 1 00. Alternatively, the cunent velocity 
program module 188 further includes a program to calculate the Doppler shift from the 
bearing angle ((p) and a radial cunent velocity, V„ measured relative to a receiver 106 in 
a backscatter configuration. 

The foregoing description of the prefened embodiments of the invention have 
20 been presented for the purposes of illustration and description. It is not intended to be 
exhaustive or to limit the invention to the precise form disclosed, and obviously many 
modifications and variations are possible in light of the above teaching. The 
embodiments were chosen and described in order to best explain the principles of the 
invention and its practical application to thereby enable others skilled in the art to best 
25 utilize the invention in various embodiments and with various modifications as are suited 
to the particular use contemplated. For example, more than two transmitters can operate 
with a single receiver. Likewise, there are situations where multiple receivers operating 
with a single transmitter may have advantages, as in military operations where greater 
radio silence is prefened. All of these possibilities constitute alternative embodiments of 
30 the present invention. Thus, it is intended that the scope of the invention be defined by 
the claims appended hereto and not limited to the precise embodiments described herein. 
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